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Abstract. We introduce a new microcanonical dynamics for a large class of Ising 
systems isolated or maintained out of equilibrium by contact with thermostats at 
different temperatures. Such a dynamics is very general and can be used in a wide 
range of situations, including disordered and topologically inhomogenous systems. 
Focusing on the two-dimensional ferromagnetic case, we show that the equilibrium 
temperature is naturally defined, and it can be consistently extended as a local 
temperature when far from equilibrium. This holds for homogeneous as well as 
for disordered systems. In particular, we will consider a system characterized by 
ferromagnetic random couplings Jij € [1 — e, 1 + e]. We show that the dynamics 
relaxes to steady states, and that heat transport can be described on the average 
by means of a Fourier equation. The presence of disorder reduces the conductivity, 
the effect being especially appreciable for low temperatures. We finally discuss a 
possible singular behaviour arising for small disorder, i.e. in the limit e — > 0. 
keywords: Transport processes (Theory), Heat conduction, Disordered systems 
(Theory) 



1. Introduction 

Microcanonical dynamics seems to be a proper mode to describe an isolated system, 
or its isolated bulk, without any assumption on the equilibrium state between the 
system and the surrounding. A typical context where such a dynamics occurs is in 
the study of transport properties in continuous or discrete models (see [1] for general 
references). For instance, for a system in thermal contact with two heat reservoirs 
at different temperatures, energy may be exchanged only through the thermostats, 
and the internal dynamics requires an energy preserving evolution to be implemented 
in the microcanonical framework; analogous considerations hold for matter transport 
and similar items [2j [3l [4]. 

Presently, we are interested in heat transport in discrete models, namely spin 
systems. As it is well known, for such systems there exist several rules to simulate 
microcanonical evolution, e.g. the Q2R rule introduced by Vichniac, the Creutz rule 
and, more recently, the KQ rule [5J El [HI [H] ■ All of them present some weak aspects 
which require the introduction of some constraints concerning geometrical support, 
boundary conditions (i.e. the allowed range of imposed temperatures), connectivity, 
homogeneity, etc. Our purpose is to introduce a microcanonical dynamics sufficiently 
flexible to be adapted to all these cases. The recent interest in systems displaying 
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quenched disorder and inhomogeneous structure, provides indeed a strong motivation 
for the development of such a general dynamics [TU1 EJ HH H31 03] ■ 

As it is well known, Q2R and Creutz dynamics differ in the fact that the latter 
supplies each site in the lattice with a possible bounded reserve of "kinetic" energy, in 
such a way that the spin flip is allowed not only when the magnetic energy is preserved 
(as in Q2R), but also when the energy excess or defect may be assigned to or extracted 
from the reserve, so that the total energy is preserved. This additional possibility, 
however, does not solve the dynamical freezing at low energy densities (occurring of 
course also in Q2R), since in such conditions one easily falls into short periodical cycles, 
inhibiting the transport of heat and the exploration of the whole set of configurations 
at constant energy. In other terms, both Creutz and Q2R dynamics have an ergodicity 
breakdown at low energy, loosing the independence of initial conditions which is 
necessary to ensure a statistical effectiveness to the dynamical description. Moreover, 
one should introduce ad hoc modifications when the connectivity is not homogeneous. 
As to the KQ dynamics, introduced and discussed in [HIS], even if it possesses the 
wanted ergodicity, it is not a satisfying solution to the general dynamical problem, 
since it is specifically meant for regular lattices. 

The new dynamics we are going to define actually modifies the Creutz scheme, 
only maintaining the distinction between magnetic and kinetic energies. As we shall 
see, not only it displays effective ergodicity for all practical purposes, but it can also 
be easily adapted to a wide class of spin systems without any substantial modification, 
exhibiting therefore the desired versatility. In particular, it can be implemented on 
arbitrary topologies, even non regular ones. Another aspect of this flexibility is that 
the extension to inhomogeneous systems, where the spin coupling is not a constant, 
is quite natural. Hence, a great variety of disordered models become easily accessible 
in principle; in the following we study an example of quenched disorder. 

It is worth underlining that, in all the cases considered here, the magnetic 
and kinetic energies result to be non correlated: this allows a natural definition of 
temperature at equilibrium, which, in a very natural way, depends only on the average 
kinetic energy. This approach can be consistently extended, when far from equilibrium, 
to a definition of local temperature. Now, dealing with heat transport, this implies an 
important improvement in the description of steady states, since not only diffusivity 
but also conductivity become easily computable. In particular, in studying a model 
of quenched disorder, the interesting point is to explore how inhomogeneity affects 
the thermal and transport properties. For instance, we verify that the conductivity 
becomes smaller as the disorder grows. The emergence of non analytical behaviour in 
the limit of vanishing disorder is also observed. 

In the following, we first describe in details how our microcanonical dynamics 
works (Section [2]) and we show that it is able to relax the system to well defined 
steady states (Section [3]), even in the presence of disorder (Section Hj), verifying that 
a temperature as an intensive quantity independent of the actual values of couplings 
can still be defined. Then, we focus on the transport properties (Section [5]), and on 
the effect of disorder on the conductivity (Section [6]) . The last section is devoted to 
our conclusions and final remarks ( Section [7]) • 

2. Microcanonical Dynamics 

Let us consider an Ising model defined on a generic network. In particular, Sj = ±1 
denotes the i-site spin; on each link i ~ j we introduce a local exchange interaction 
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Jij, so that the magnetic energy of a link is JijSiSj. Moreover, on the link i ~ j we 
define a local kinetic energy E t j > 0. We remark that while standard Creutz scheme 
specifies a bounded amount of "kinetic" energy in each site, we introduce an energy 
defined on links which is non-negative but, in principle, unbounded. The dynamical 
rule proceeds as follows: 

(i) Start from a distribution of energies ; 

(ii) choose randomly a link i ~ j; 

(iii) extract one over the possible four spin-configurations for the couple of sites i,j, 
and evaluate the magnetic energy variation AE m induced by the move, where 
AE m allows for energy variations occurred on the link i ~ j as well as on those 
pertaining to links adjacent to sites i or j; 

(iv) if AE m < accept the move and increase the link energy Eij of AE m . When 
AE m > 0, accept the move and decrease the link energy of AE m only if 
Eij > AE m . 

Three main points have to be remarked: 

• Every link is defined by its extremes not only in rectangular lattices but in all 
graphs, even in the case of variable connectivity. This ensures a great flexibility 
with respect to the underlying geometry. 

• The interchange between sites and links allows the energy to travel throughout the 
graph, even at low energy density, without the freezing of the standard Creutz- 
Q2R rules or the geometrical restrictions of the KQ dynamics. 

• The constraint of uniform couplings Jij = J = const, typical of previous 
microcanonical dynamics, can also be relaxed as our dynamics allows to use local 
(i.e. non uniform) spin interactions. 

Clearly, the above dynamics conserves the total energy given by the following 
Hamiltonian function 



where the sum runs over all the links of the network. Notice that energy is locally 
conserved, namely long-range energy transfers are forbidden and hence we can define 
our dynamics as microcanonical, where the kinetic energy energy E^ works as an 
additional degree of freedom. Expression JT]) is very general, representing a conserved 
energy for a model defined on a generic network. Moreover, according to the choice 
of the couplings Jij, a wide class of models (disordered interactions, vacancies, spin 
glasses etc.) can be recovered. 

Actually, for comparisons and checks, we consider first the homogeneous system 
Jij = J (with J = 1 in simulations) on the standard rectangular lattice of Lx 
columns and Ly rows, both in the toroidal geometry (periodical boundaries) and 
in the cylindrical geometry with two opposite boundaries. Such boundaries can be 
"free" (i.e. three neighbours for each site) or "open" (i.e. in contact with a reservoir). 
Heat conduction in the latter case was studied e.g. in [T^l HB] using Q2R-Creutz 
dynamics, in [9] using KQ dynamics and in [17] using a layer heat-bath. 

Afterward, we shall study the transport properties in a cylindrical system where 
the couplings {Jij} are chosen randomly in the interval [1 — e, 1 + e], i.e. a disordered 
ferromagnet. 




(1) 
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3. Temperature and Dynamical checks 

At equilibrium, i.e. in a closed system where no energy is injected or subtracted, 
every move is reversible, and the probability P(S — > S') to go from the configuration 
S to the configuration S' equals P(S' — > S). Detailed balance is satisfied, and the 
equilibrium probability distribution is constant over all configurations. 

If all the configurations with the same energy are visited in the evolution, the 
system is ergodic and its equilibrium properties can also be described by Boltzmann 
statistics. In this perspective, we check the agreement between the numerical 
stationary results and the Boltzmann distribution, focusing first on the homogeneous 
system = J defined on a two-dimensional torus; a positive response to the check 
has to be seen as an indication of effective ergodicity. 

A first step in this direction consists in testing the decorrelation between the 
kinetic energy 25y and the magnetic energy JsiSj of a link. Indeed, we verify that for 
each link 

(EijJsiSj) = (Eij)(JsiSj) , (2) 

where (■) denotes a time average along our dynamics. Kinetic and magnetic part of 
the Hamiltonian can therefore be treated separately. 

As a second point, we observe that the link energy satisfies the Boltzmann 
distribution exp(— /?.Ey), where the fitted constant (3 = 1/T is naturally the inverse 
temperature of the system (Fig.[TJ left-panel). In general, for the homogeneous system 
we can write E%j = AJnij + £y- , where n%j is an integer and £jj is a constant depending 
on the initial conditions; since £y does not have any influence, hereafter we set £y = 0. 
Notice that the factor 4 multiplying Jn^ is due to the fact that the minimal amount 
of exchanged energy is 4J {Eij derives from AE m and, apart from possible borders, 
in rectangular lattices the coordination number is even). 

The equilibrium temperature has been verified to be independent from the initial 
distribution of energies and from the particular link considered, as expected from an 
intensive quantity. Averages on the links are possibly used only to accelerate the 
numerical convergence. 

The simplest way to evaluate the temperature from the above distribution is to 
relate (3 with the average link-energy, obtaining 

V°° 4 7r?p~ 4J ' 3 " 4 J 

\ ' — e -4J/3n e 4/3J _ I \°) 

which may be easily inverted to get (3: 

=h l ^{ l+ w))- (4) 

In a subsequent check, we verify that T = l/j3 given by formula (0]) is the 
temperature even for the magnetic variables of the system. Indeed, Fig. [2] evidences 
that the average magnetic energy (E m ) and its fluctuations ((E m ) 2 ) — (E m ) 2 are 
the same in the microcanonical system at temperature T given by Eq. [4] and in a 
two-dimensional Ising model in equilibrium at temperature T. 

Equivalent results are also obtained introducing thermostats at the boundaries of 
the microcanonical system (FigJ^]). The coupling between system and thermostats is 
imposed by extracting the link energies at the borders every time-step according to the 
probability distribution e~ 4/3J ™. We remark that in presence of such thermostats the 



Energy Transport in an Ising Disordered 



Model 



5 




Figure 1. Probability distribution of the kinetic energy per link for an Ising 
system defined on a square torus of size Lx = Ly = 48 with constant coupling 
strength J (left panel) and with random J and disorder e = 0.2 (right panel). 
Results found for different initializations Eg are reported as shown by the legend. 
Notice that in the right panel the continuous energy distribution p(E, Eg), J?€l 
is approximated by the histogram p(n, Eg), n E N. 



whole system equilibrates at temperature 1/(3, independently of the initial conditions, 
proving the reliability of such a simple model for thermostats. 

Energy transport is naturally introduced by imposing different temperatures at 
the boundaries. Then, whenever a stationary state is reached, Equation Q provides 
a simple definition of local temperature. Clearly, along the same column all energies 
and temperatures are expected to be equal, as they are in fact. The consistency of this 
definition of local non-equilibrium temperature has been further verified by showing 
that, at the same T, the local magnetic energy for a link is the same at equilibrium 
or in the presence of transport (Fig. \2§ ■ 

In conclusion, the homogeneous system follows a quasi equilibrium picture 
consistent with the results given by the KQ dynamics [9]. The remarkable difference 
is that, in the present case, the existence of a local kinetic energy allows a natural 
definition of temperature, without any reference to the equilibrium magnetic model. 

4. From Homogeneous to Disorderd Lattices 

After testing our dynamics for the Ising model on a regular lattice with homogeneous 
coupling, we exploit the fact that, as anticipated, this dynamics is naturally defined 
even in more complex situations, e.g. when the underlying topology or the interaction 
couplings Jij are not homogeneous. Here in particular we focus on the case of a 2- 
dimensional rectangular lattice where the coupling Jjj are randomly chosen from a 
uniform distribution in an interval, i.e. Jij € [1 — e, 1 + e], < e < 1; hence, e is an 
index of disorder, being the average magnetic interaction fixed and equal to 1. 

First of all, we study the differences at equilibrium between homogeneous and 
disordered systems. To this purpose, using the standard Metropolis dynamics, we 
evaluate the critical temperature Tc by means of Binder cumulant techniques [18J. 
Figure (J3J) shows the value of the critical temperature as a function of the disorder e. 
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Figure 2. Average magnetic energy per spin (E m ) (main figure) and its 
fluctuations ((E m ) 2 ) — (E m ) 2 (inset) as a function of the temperature T for 
Ising systems of size TV = 2304 defined on a torus and analyzed with the 
Metropolis canonical dynamics (continuous line) and with the microcanonical 
dynamics introduced here (dark •); the case of an Ising system defined on a 
cylinder endowed with two thermostats set at the same temperature (bright •) or 
at different temperatures (x) is also considered. The consistency among all cases 
is clear. 



We note that Tc is slightly lowered by the presence of disorder. 

Let us now consider the new dynamics for the disordered system. With respect 
to the homogeneous case Jij = 1, the main difference is that the link energy is not 
of the form Eij = Amj with mj integer. Indeed, the kinetic energy on the link i ~ j 
can be written as Eij = Y^k~h Jkhnkh, where Eij > 0, nkh are integers and the 
sum X)fc~h runs over i ~ j and its six neighbouring links. Hence, in the sum there 
are seven different real numbers Jhk randomly chosen in [1 — e, 1 + e] and, if all the 
combinations of integers n^h are realized in the dynamical evolution, the values of E^ 
can be considered densely distributed along the positive real axis (see Appendix). 

Therefore, under the hypothesis of ergodicity we expect that the probability 
density for link energies is exp(— @E), with E a real number. In Figure (|TJ) we 
check indeed that in a random system the link energies are distributed according 
to a continuous exponential function. Also in this case, as in ([3]), the temperature 
can be related to the average kinetic energy, in particular for a continuous variable we 
have 



(E) = 



J °° z e @ z dz 



1 

J 



(5) 
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Once again, the constant 1//3 is the temperature T and the link (kinetic) energy equals 
the system temperature. We check that, although the system is inhomogeneous, the 
temperature is the same for all the sites, as expected at equilibrium, and also that (E) 
is a consistent definition of temperature even for the magnetic system, as analogously 
verified for the homogeneous system (see Fig. [2]). Therefore, the temperature is the 
parameter governing the equilibration of the system. In particular, at equilibrium 
different parts of an inhomogeneous system may display different behaviours but T is 
constant in the whole sample. For example, if a part of the system is characterized 
by constant coupling J,j = J while in the rest Jij is chosen randomly, the local 
kinetic energy is described by discrete and continuous distributions in the different 
regions respectively. However, the temperature is the same in both. Analogously in 
a free cylinder (isolated borders) with constant coupling J, we have that the border 
sites exchange energy with an odd number of neighbours. The minimal amount of 
exchanged energy is therefore 2 J instead of 4 J, giving a different discrete distribution 
for sites in the bulk and at the border. However, we have numerically verified that 
the temperature does not change. 

Thus, in general, for systems with random-incommensurate couplings, Equation 
(|5j) provides a good definition of temperature. This definition can naturally be 
extended to the case when different temperatures are imposed at the boundary 
thermostats, i.e. in the presence of gradients (T% — T± > 0) and heat transport, 
as described in the next Section. 
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5. Transport, Conductivity, Diffusivity 

As previously explained, our system satisfies the Boltzmann distribution exp (-E/T), 
so that in order to fix a temperature T on a given set of links, it is sufficient to 
extract the relevant link energies according to the distribution exp (-E/T), and this is 
sufficient to establish a consistent coupling between the system and the thermostats. 
In the following we consider a cylindrical Ising system with boundaries coupled to 
two thermostats at temperatures 7\ and T 2 , respectively, and we study the transport 
properties. First of all, we verify that in the steady state the kinetic energy (-Ey) is still 
well defined on each link. However, in such a non-equilibrium situation, it varies from 
link to link; hence, hereafter we use (-Ey) as a definition of local temperature Ty. We 
also notice that once the distance from the thermostats is fixed, its spatial fluctuations 
increases with e. Hence, the dependence displayed by Ty on the magnetic coupling 
Jij is very small for e not too large. Now, in the presence of temperature fluctuations 
also transport is expected to display non trivial spatial patterns arising from the 
interplay between non-equilibrium and disorder. As a first step in studying such a 
complex phenomenology we investigate how disorder affects the average transport 
properties. In particular we introduce the column temperature Tj and magnetic 
energy EJ 1 obtained by spatially averaging over all i's Ty and (E^), respectively. 
The emerging cylindrical symmetry makes j the only gradient variable, and in this 
situation the simplest model for heat transport is provided by the one-dimensional 
Fourier equation, where, as usual, the column index j will be treated as a continuous 
variable a;: 

dE tot Q f dT \ 

"«r ■ ( 6 ) 



dt dx \ dx ) 

Notice that E tot is the specific total energy in column j, namely Ej 0t = Tj + EJ 1 , and 
that k is the conductivity which is expected to depend only on the local temperature 
and on the system parameters (J and e in this case). Indeed, in the following we verify 
that energy transport can be studied by means of equation and we also evaluate 
k as a function of T, evidencing the effect of disorder. 

In the stationary state, when dE tot /dt = 0, we call Q the total amount of 
heat crossing any given surface, so that dQ/dt represents the heat flux which can be 
measured, for instance, from the energy absorbed per unit of time by the thermostat. 
As shown in [8], we have 

clearly k represents the proportionality constant between the temperature gradient 
and the flowing heat energy. Similarly, the diffusivity of the magnetic energy can be 
defined as: 

dQ/dt 

dE m /dx V ' 

It is worth underlining that while D can be obtained directly from the magnetic energy 
E m (see Eq. [8]) for the conductivity k we need a measure of the temperature T (see 
Eq. [?])■ For this reason KQ dynamics, dealing exclusively with the magnetic energy, 
allows for the computation of diffusivity, while here we are able to easily estimate both 
quantities. 
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Figure 4. Conductivity versus temperature for an ordered (constant J) Ising 
system. Different sizes L and different temperature gaps AT are considered, as 
shown in the legend. 



Introducing the magnetic specific heat C = dE m /dT, a comparison of Eq. [7] and 
Eq. [5] provides the simple relation 



Notice that C is an equilibrium quantity well known for the Ising model. 

We have measured conductivity for the ordered system with e = (Fig. [J} as well 
as for the disordered system with fixed e ^ (Fig. [5]). In both cases we found a very 
good collapse of data pertaining to different temperatures intervals and sizes. 

This provides a strong confirmation of the picture above based on a Fourier 
transport equation with a temperature (energy) dependent conductivity (diffusivity) . 
In particular, we notice that the curve for the conductivity exhibits a maximum at a 
temperature T slightly larger than the critical temperature Tc which corresponds to a 
peak in the specific heat and to a minimum for the diffusivity (see Fig. [6]). As observed 
also in [9], we expect that in the thermodynamic limit the diffusivity goes to zero 
logarithmically, consistently with the logarithmic divergence of the specific heat, so 
that the conductivity is not expected to diverge at the critical point To- Interestingly, 
from both Fig. [J] and Fig. [5J we notice that as the temperature is increased, the rate 
of growth (for T < T) and of decrease (for T > T) displayed by the conductivity are 
sensitively different. This point can be understood by looking at Fig. [6] and recalling 
Eq. [9] at small temperatures the diffusivity remains finite, while the specific heat is 
close to zero. 
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Figure 5. Conductivity versus temperature for a disordered (e = 0.2) Ising 
system. Different sizes L and different temperature gaps AT are considered, as 
shown in the legend. 



6. Effects of Disorder 



Let us now deepen our knowledge of how the conductivity is affected by disorder 
e. In Fig. [7] we show data for conductivity vs. temperature for different degrees of 
disorder. Remarkably, the figure highlights a clear drift: as the disorder e increases the 
conductivity k gets lower for any value of the local temperature. Hence, for a given 
system characterized by the set of parameters (Lx, Ly , Ti, T 2 ) the heat flow Q(e) is 
monotonically decreasing with e, or, otherwise stated, the presence of disorder lowers 
the transport efficiency. We now focus on the behaviour for e — * 0. The transition 
from the ordered to the disordered case is in fact characterized by different expressions 
of the kinetic energy as a function of the temperature (see Eqs. HI [S]). This could be 
a signal of a non trivial behaviour also for different physical quantities. A plausible 
ansatz for a singularity at sufficiently small e is provided by the expansion: 

Q(e)~Q(0)-Ae a , (10) 

where A > is expected to depend in general on the system parameters (size, gradient 
Ta — Ti, etc) and a is a possibly universal exponent. In Fig. [5] we show a log-log scale 
plot for Q(Q) — Q{e) which confirms the scaling law of Eq. ITTJ1 over a wide range of 
disorder and for different choices of temperature intervals. 

By means of a fitting procedure it is possible to get estimates for A and a. 
Interestingly, A increases with the temperature gap T 2 — Tl , and its value is significantly 
smaller when both temperatures are larger than the critical one, namely T 2 , Tl > Tc- 
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Figure 6. Specific heat C (continuous line) and diffusivity D (dashed line) as a 
function of the temperature T, for an ordered Ising system (e = 0) of size L = 16 
and thermostats at temperatures T\ = 0.01, T2 = 5.0. 



Hence, as expected, the effect of disorder is especially appreciable for low temperatures. 
We also notice that, within the error, the agreement among the exponents a is rather 
good. 

Finally, we remark that when e <C 1 and both T\ and T2 are smaller than Tc the 
decorrelation time are long and the measures of Q(e) and a, the relevant exponent, 
get awkward. Indeed, the "critical slowing down" observed for e — > is neatly distinct 
from the behaviour at J = 1 (i.e. e = 0). This provides a hint on a possible non 
analiticity as e — ► 0. 

7. Conclusions 

We have introduced a microcanonical dynamics based on the link reserves of kinetic 
energy, showing that it recovers the results on regular lattices previously obtained by 
other dynamics and that, moreover, it can be properly applied to disordered systems 
and it lets a natural definition of temperature, also when far from equilibrium. This 
allows an easy access to such quantities as conductivity and specific heat, otherwise 
difficult to define operationally. We measured the average conductivity for the ordered 
and the disordered systems showing that the presence of disorder reveals in the 
reduction of the transport efficiency. 

The disordered lattice exhibits remarkable additional features: there are 
indications for a singular behaviour of the system in the limit e — > of the disorder 
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Figure 7. Conductivity versus temperature for a disordered Ising system of size 
L = 16 and fixed temperature gap T\ = 0.01, T2 = 5.0. Different values of disorder 
e are considered, as shown in the legend. 

parameter. The reason of such a singularity may lie in the numerical strain in filling the 
reals by combinations with integer coefficients if the couplings are very close to each 
other. This slowing down should be the dynamical counterpart of the discontinuity 
of the temperature (see expressions ([3]) and J5])) when the 2e- width of the disorder 
interval tends to zero. Moreover, it would be interesting to highlight the nature of 
the exponent a, which describes how the heat flow decreases with e. In particular, it 
should be clarified whether a depends on the system parameters (Lx, Ly, T2) or, 
if not, whether it depends on the particular dynamics chosen. 

Recalling from past experiments that the sole presence of gradients proved influent 
on the fluctuations of physical and geometrical observables [9] , on the basis of present 
results it seems likely that there are other problems deserving future investigations, 
e.g. the fine interplay between disorder and non-equilibrium in small scale regions, 
where local heat flows and temperature fluctuations interact presumably in non trivial 
ways. In other words, for the disordered system we have verified the effectiveness of 
the Fourier approach by considering averaged quantities, while the validity of such a 
description at the microscopic scale is still an open problem. 

Finally, we underline that the dynamics introduced here allows possible extensions 
not only to varied types of disorder (vacancies, negative couplings, dynamical 
disorder...) or to spin systems lying on different geometrical structures, but also 
to systems where the gradient maintaining the condition of non-equilibrium is due, 
for instance, to the contact with reservoirs at different densities [T5]. In general, it is 
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conceivable to model some peculiar dynamics where "conservation" refers to the sum 
of two distinct populations. 

Appendix A. 

Let 9i,...,9 m , with m > 2, positive and rationally independent reals. Let X > 0. 
Then, for arbitrarily small 6, there exists a linear combination with integer coefficients 

m 

C = N k 6 k 
fc=i 

such that <| X - C |< 5. 

Proof : let Ao and Mq be the minimal and the maximal positive reals obtained as 
integer combinations X)fcLi u fc^fcj where u k may be 1,0,-1. Then define ao as the 
greatest positive real < ao < Ao of the form 

a = M - fc A , 

with integer fco- Such a number ao is an integer combination of the {9k} and it does 
not coincide with Ao because of their rational independence. We define now 

Ai = min{ao, Ao - ao}, 
observing that necessarily Ai < Aq/2. 
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Now, we iterate the procedure denning 

a n = Mo — k n \ n , 

and 

K+i = min{a„, A„ - a n }. 

We obtain a positive sequence A„ converging to 0, because A n +i < A n /2. Each term 
in this sequence is an integer combination of the {9k}- Now, we can choose N in such 
a way that Ajv < S. Since Ajv is finite, there exists an integer p such that 

p X N < X < (p+ 1) Ajy. 

We get this way a neighborhood of X smaller than 5, whose extremes are integer 
combinations of the {8k}, QED. 

Note 1: the way X has been approximated does not necessarily correspond to an 
optimal strategy; 

Note 2: with reference to our problem, where the {9k} are the seven couplings 
chosen for each link in an interval [1 — e, 1 + e], and the integer coefficients are the 
iterated additions and subtractions determined by the wandering spin flips, it seems 
that, especially for great X'b, the probability of the right approximating combination is 
extremely low. However, in the integrals of formula (|5j). the real variable is weighted in 
turn by the Boltzmann factor, and this explains why also in our finite-time simulations 
it is possible to to fill densely the relevant part of the real axis, ensuring the correctness 
of the procedure; 

Note 3: when the 2e-width of the interval shrinks (little disorder), this contributes 
to make difficult the approximation: the reduction in numerical efficiency in these 
conditions may be observed indeed, as noticed in Section |5]). and it could also be 
related to the possible singularity of the limit e — > 0. 
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